Spin polarization dependence of quasiparticle properties in graphene 
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We address spin polarization dependence of graphene's Fermi liquid properties quantitatively using 
a microscopic Random Phase Approximation theory in an interacting spin-polarized Dirac electron 
system. We show an enhancement of the minority-spin many-body velocity renormalization at fully 
spin polarization due to reduction in the electron density and consequently increase in the interaction 
between electrons near the Fermi surface. We also show that the spin dependence of the Fermi 
velocity in the chiral Fermi systems is different than that in a conventional two-dimensional electron 
liquid. In addition, we show that the ratio of the majority- to minority-spin lifetime is smaller than 
unity and related directly to the polarization and electron energy. The spin-polarization dependence 
of the carrier Fermi velocity is of significance in various spintronic applications. 
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I. INTRODUCTION 

Graphene is a two-dimensional crystal of carbon atoms 
has been recently discovered^ This stable crystal has 
attracted considerable attentiorP because of its unusual 
effective many-body properties"^ that follow from chiral 
band states and because of potential applications. The 
low energy quasiparticle excitations energy in graphene 
are linearly dispersing, described by Dirac cones at the 
edges of the first Brillouin zone. 

Stable non-reactive graphene layers on top of ferro- 
magnetic materials 9 might be used as sources of spin- 
polarized electrons. Electron sources are used in all do- 
mains ranging from technical devices like cathode-ray 
tubes to large scale scientific experiments like electron 
accelerators. This is of great interest for studies of mag- 
netic systems in condensed matter physics, including the 
field of spintronics. 

Graphene's spin-transport properties are expected to 
be particularly interesting, with predictions for extremely 
long coherence times and intrinsic spin-polarized states 
at zero fielcPn Spin-polarized electron emission from the 
graphene/Ni system before and after exposure to oxygen 
has been recently studied^ and the spin polarization of 
secondary electrons obtained from this system upon pho- 
toemission and suggested to use such passivated Ni sur- 
face as a source of spin-polarized electrons stable against 
adsorption of reactive gases. The resolve spin transport 
from conductance features that are caused by quantum 
interference has already been measured^. These features 
split visibly in an in-plane magnetic field, similar to Zee- 
man splitting in atomic and quantum-dot systems. As a 



result, spin- up and spin-down conductance contributions 
at finite field arc offset in gate voltage, leading to Zeeman 
splitting of interference features in a gate-voltage^. 

Many electronic and optical properties of graphene 
could be explained within a single-particle picture in 
which electron-electron interactions are completely ne- 
glected. The discovery of the fractional quantum Hall 
effect in graphene^ represents an important hallmark in 
this context. By now there is a large body of experimen- 
tal worlPEEHU showing the relevance of electron-electron 
interactions in a number of key properties of graphene 
samples of sufficiently high quality. 

Conventional two-dimensional electron gas (2DEG), on 
the other hand, has been a fertile source of surprising new 
physics for more than four decades. Although graphene 
was only isolated for the first time in 2004 and it is still at 
an early stage, it is already cleai^ that the strong field 
properties of Dirac electrons in graphene are different 
from and as rich as those of a semiconductor hetero junc- 
tion 2DEG. The Fermi liquid phenomenology of Dirac 
electrons in graphene^' and conventional 2DEGp^ have 
the same structure, since both systems are isotropic and 
have a single circular Fermi surface. The strength of in- 
teraction effects in a conventional 2DEG increases with 
decreasing carrier density. At low densities, the quasipar- 
ticle weight Z is small, the velocity is suppressed^, the 
charge compressibility changes sign from positive to neg- 
ative, and the spin-susceptibility is strongly enhanced^. 
These effects emerge from an interplay between exchange 
interactions and quantum fluctuations of charge and spin 
in the 2DEG. 

In addition, effective mass or the effective Fermi ve- 
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locity is an important concept in Landau's Fermi liq- 
uid theory since it provides a direct measure of the 
many-body interactions in the electron system. In the 
highly interacting, dilute, paramagnetic regime in the 
2DEG the effective Fermi velocity, which is defined by 
the effective mass as v* — hkp/m*, is significantly di- 
minished compared to its band value and tends to de- 
crease with increasing rj^^l, the so-called Wigner-Seitz 
radius. Recent measurements of the effective mass for 
two-dimension electrons confined to AlAs quantum wells 
revealed that, when the 2DES is fully valley- and spin- 
polarized, the effective mass is suppressed down to values 
near or even slightly below the band mas d 19 * 2 ^ ^. A so- 
phisticate theoretical calculation has been showrP^l that 
in an interacting, fully spin-polarized 2DES the absence 
(freezing out) of spin fluctuations reduces the effective 
mass below its band value, in agreement with experi- 
mental data. Furthermore, the spin-up and spin-down 
effective masses from magnetotransport measurements at 
different temperatures for a 2DEG and the effective hole 
mass measurements through analyzing the temperature 
dependence of Shubnikov-de Haas oscillations in dilute 
2D hole systems have been recently reported^. 

In the Dirac electrons in graphene, it was showrP^H 
that interaction effects also become noticeable with de- 
creasing density that the quasiparticle weight Z tends to 
larger values, that the velocity is enhanced rather than 
suppressed, and that the influence of interactions on the 
compressibility and the spin-susceptibility changes sign. 
These qualitative differences are due to exchange interac- 
tions between electrons near the Fermi surface and elec- 
trons in the negative energy sea and to interband contri- 
butions to Dirac electrons from charge and spin fluctua- 
tions. 

Our aim in this work is to study the spin polariza- 
tion dependence of quasiparticle properties in graphene 
particulary the renormalized velocity and inelastic scat- 
tering lifetime of quasiparticles within the leading-order 
single-loop self-energy expansion. Our theory for spin 
polarization dependence of quasiparticle velocity renor- 
malization in interacting Dirac electron systems is mo- 
tivated not only by fundamental many-body considera- 
tions, but also by application to improve high-speed op- 
eration in the spintronic devices^ and potential future 
experiments. By chemical doping in graphene, spin po- 
larization effects are predicted for some adsorption con- 
figurations^!. Remarkably, the studies of spin polariza- 



tion dependence of quasiparticle properties should help to 
understand spin valve physics and recent measurements 
of electronic spin transport in garphene^, and the possi- 
bility of magnetism in graphene induced by single carbon 
atom defectsP. 

The paper is organized as the following. In Sec. II 
we introduce the formalism that will be used in calcu- 
lating spin polarization quasiparticle properties which 
includes the many-body effects by suing RPA. In Sec. 
Ill we present our analytical and numerical results for 
the self-energy and renormalized Fermi velocity in doped 
graphene sheets. Sec. IV contains discussions and con- 
clusions. 

II. METHOD AND THEORY 

We consider the long-range Coulomb electron-electron 
interaction. We left out the intervalley scattering and use 
the two component Dirac Fermion model. Accordingly, 
the total interacting Hamiltonian in a continuum model 
at K + point is expressed aP^ 

H = -iHvJ2*-Vi + l^2V(r i -rj), (1) 

i i^j 

where a are Pauli matrices and v — 3ta/2h ~ I0 6 
m/s is the Fermi velocity with a ~ 1.42 A is the carbon- 
carbon distance in honeycomb lattice. Here p 4 = —ihVi 
is the canonical momentum of the i— th electron and v q = 
27re 2 /eq is the Fourier transform of the bare Coulomb 
interaction where e is an average dielectric constant of the 
surrounding medium. The coupling constant in graphene 
or graphene's fine-structure constant is a ee — e 2 /ehv. 
The coupling constant in graphene depends only on the 
substrate dielectric constant while in the conventional 
2D electron systems is density dependent. The typical 
value of dimensionless coupling constant is 0.25 or 0.5 for 
graphene supported on a substrate such a SiC or Si02- 

As it is clearly seen from the first term of Eq. ([I]) , the 
spectrum is unbounded from below and it implies that 
the Hamiltonian has to be accompanied by an ultraviolet 
cut-off which is defined k c and it should be assigned a 
value corresponding to the wavevector range over which 
the continuum model Eq. ([I]) describes graphene. For 
definiteness we take k c to be such that 7rfc 2 = 2(27r) 2 / Aq, 
where Aq = is the area of the unit cell in the 

honeycomb lattice. With this choice, the energy hvk c = 7 
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eV and 
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The continuum model is useful when fc c 3> kp , i. e. when 
A ^> 1. Note that, for instance, electron densities n = 
0.36 x 10 12 and 0.36 x 10 14 cm' 2 correspond to A = 100 
and 10, respectively. 

The spin-polarization dependence of dynamical po- 
larizability tensor in terms of one-body noninteracting 
Green's function is written as^ 



Tt[i 1o G^ (k + q, u + n, M)*7oGi 0) (k, w, m)] (3) 

where tr refers to the spin-direction, j" or J,. After im- 
plementing (k, oj, jLt) in Eq.|3| and calculating the in- 
tegral, the results end up to the follow expressiorP 
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where g v = 2 is valley degeneracy. fi a is the spin depen- 
dence chemical potential, B = q 2 /(8 ir^/Vt 2 + v 2 q 2 ) and 

c= (2/i CT + %n)/{vq). 

The technical calculation^! on which our conclusions 
are based is an evaluation of the spin-polarization depen- 
dence electron self-energy SJ(k, u) of the Dirac fermion 
near the quasiparticle-pole. £f(k, lo) describes the in- 
teraction of a single Dirac electron with spin a near the 
2D Fermi surface with all states inside the Fermi sea, 
and with virtual particle-hole and collective excitations 
of the entire Fermi sea. A direct expansion of electron 
self-energy in powers of the Coulomb interaction is never 
possible in a 2D electron liquid because of the long-range 
of the Coulomb interaction. Our results for the Dirac 
electron gas are based on the random phase approxima- 
tion (RPA) in which the self-energy is expanded to the 
first order in the dynamically screened Coulomb interac- 
tion (setting h = 1): 
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G°?(k + q,iu; n +ffi m ),(5) 



where s = + for electron-doped systems and s = — for 
hole-doped systems, £ is the spin polarization parameter, 



C = 1% - ni\/n, f3 = l/(k B T) and e(q,iCl m ,C) is the 
RPA dielectric function. n„ is the spin polarized density 
and n is the total density of system. The RPA dielectric 
function is given by 



£ (q,»o,C) = i-«Jxf ) (q,^,C) + xr(q^^,C)] ■ (6) 

In the Dirac 2D electron system, dielectric-function 
contributions from intraband and interband excitations 
are subtly interrelated. The two contributions must be 
included on an equal footing in order to describe the 
Dirac fermion physics correctly. For example one key 
property, that the static dielectric function is indepen- 
dent of q at small q, requires intraband and interband 
contributions to be summed. 

For definiteness, we limit our discussion to an electron- 
doped system with spin polarization dependence of chem- 
ical potential \i° . 

In Eq. Q u n = (2n + 1)tt//3 is a fermionic Matsubara 
frequency, the sum runs over all the bosonic Matsubara 
frequencies fi m = 2mw//3. The factor in square brackets 
in Eq. |5]), which depends on the angle #k,k+q between 
k and k + q, captures the dependence of Coulomb scat- 
tering on the relative chirality ss' of the interacting elec- 
trons. The Green's function G° s a {k 7 iuj) = ~ fj(k)] 
describes the free propagation of states with wavevector 
k, Dirac energy £J(k) = svk — /i CT (relative to the chem- 
ical potential) and chirality s — ±. The quasiparticle 
exciatation energy measured from the chemical poten- 
tial can be given by Dyson equation E% (k) = £s(k) + 
5ie£i ret ' <T ' ) (k, ui) evaluated at u> = El (k). After continua- 
tion from imaginary to real frequencies, iu — > u> + irj and 
using the Dyson equation, the spin dependence renor- 
malized Fermi velocity can be expressed^! in terms of the 
wavevector and frequency derivatives of the retarded self- 
energy J]^ ct ' cr ^ (k, uj) evaluated at the spin dependence 
Fermi surface which is k a F = (1 + oQY^kp, where kp is 
Fermi momentum: 
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FIG. 1: (Color online) Spin-polarization dependence renor- 
malized velocity, in the Dyson scheme, scaled by that of a 
noninteracting velocity as a function of the degree of spin po- 
larization C for cut-off value A = 100 (n = 0.36 x 10 12 cm" 2 ) 
in (a) the down-spin where < ni < 0.18 x 10 12 cm" 
(b) the up-spin where 0.18 x 10 12 < n t < 0.36 x 10 12 cm 
for different coupling constant values. The spin polarization 
dependence of the up- and down-spin velocity behaves differ- 
ently as £ increases. 



hand, the renormalized velocity is given by 

*(OSA) 
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This expression can also be obtained from the formal 
definition of v* given in the first equality in Eq. (|7|) when 



the second term in the Dyson equation, 5ie£+ et (k, w), 
is evaluated at the bare pole, u = £+(k). The OSA 
thus gives the quasiparticle velocity to the first order in 
the retarded self-energy. The renormalized velocity in 
this approximation demonstrates qualitatively the same 
behavior obtained by the Dyson equation, Eq. |7]) but its 
magnitude is larger than the one calculated within the 



Dyson scheme. 

The quasiparticle weight factor Z a evaluated at the 
spin dependence Fermi surface and given by Z~ x — 

In the up-spin case, the 



1- 9 w 3?eSV ret ' ff) (k,w) 

majority-spin, Z-\ value is a bit smaller and the down- 
spin case, the minority-spin, Z± value is bigger than the 
results of Z(Q = 0). 



III. NUMERICAL RESULTS 

Since the single-particle self-energy, the density of 
states, the dynamical screening, the Fermi momentum, 
and the Fermi energy in the chiral Dirac fermion are all 
affected by spin polarization, we expect all Fermi liq- 
uid parameters to be strongly dependent on the spin- 
polarization parameter. An important thermodynamic 
quantity is the system compressibility, which has been 
already studied by two of uJ^. 

Our results for spin-polarization dependence of the 
Dirac electron velocity, v* a /v, at fixed electron den- 
sity value in the up- and down-spin, the majority- and 
minority-spin as a function of the £ are summarized 
in Fig. 1 for different values of the dimensionless cou- 
pling constant, a ee . For the up-spin Dirac electron, 
renormalized velocity decreases with increasing spin- 
polarization degree of freedom. However, the down-spin 
electron renormalized velocity increases by increasing 
spin-polarization. These behaviors are based on the effect 
of the exchange energy in the spin channels between elec- 
trons near the Fermi surface. Since the electron density 
in the down-spin channel, = n(l — Q/2 is less than the 
electron density in the up-spin channel = n(l + £)/2, 
and decreases however increases by increasing Q, 
therefore, the exchange contribution of the down-spin is 
dominated and results in increasing the renormalized ve- 
locity in the spin-down channel and decreasing the renor- 
malized velocity in the spin up-channel. In contrast, the 
2DEG where the down-spin mass increases with spin- 
polarization first and as £ approaches near to one, it 
decreases sharply, the spin-polarized down-spin velocity 
tends to a constant when £ reaches to unitjW. It should 
be noticed that the up- and down-spin Fermi velocities 
are the same value at ( = 0. 

As it has been discussed previousljE^, graphene's 
Fermi liquid properties depend only weakly on the carrier 
density which is expressed in terms of the cut-off param- 
eter A. The trends exhibited in Fig. 1 can be understood 
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FIG. 2: (Color online) Renormalized velocity of the down-spin 
scaled by that of a noninteracting velocity, f|/i>, as a function 
of the coupling constant a ee for cut-off value A = 100 (n = 
0.36 x 10 12 cm -2 ) in (a) Dyson and (b) OSA approximations 
given by Eqs. Q and ([Sj, respectively. 



FIG. 3: (Color online) Renormalized velocity of the up-spin 
scaled by that of a noninteracting velocity, /v, as a function 
of the coupling constant a ee for cut-off value A = 100 (n = 
0.36 x 10 12 cm -2 ) in (a) Dyson and (b) OSA approximations 
given by Eqs. |7| and QSj , respectively. 



by considering the limits of small a ee and the limit of 
large q at all values of a ee . In the former limit screen- 
ing is weak except at extremely small q. The self-energy 
can be decomposed as the sum of a contribution from 
the interaction of quasiparticles at the Fermi energy, the 
residue contribution E( res ' CT \ and a contribution from in- 
teractions with quasiparticles far from the Fermi energy 
and via both exchange and virtual fluctuations, the line 
contribution £( lmo >< T ). In 9 td E J ^ cs,CT ' ) (k, uj), for example, 
the integral over q diverges logarithmically at small q 
when e(q, uj = 0) is set equal to one, i.e. when screening 
is neglected. Accordingly, screening cuts off this loga- 
rithmic divergence at a wavevector. More precisely, we 



have 
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and notice that 9 (J 3?eS_^ cs ' <T ' ) (hp, 0) = for a case that 
cjQ = —1. Because e(q, uj — 0) happens to be indepen- 
dent of q for transitions between Fermi surface points, it 
is possible to evaluate 9 W S_^ CS ' <T ' (k, uj) for the case that 
aQ 7^ — 1 analytically We find that 

\k Fl uj = 0) = — x 



y/4 - T)l In 



(10) 



where T\ a = g v a ee (^/l + ( + \/l - C)/\/l + oQ. It should 
be worthwhile mentioning that for r\ a > 2 we use an 
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equality in which —i \n(x + i\J\ — x 2 ) = arccos(a;) when 
x < 1. Importantly, similar £ dependence appears in 
the v* at small a ee . More precisely, the C dependence 
of u*(C) , both in the Dyson and OSA schemes, in the 
a ee —¥ limit, for £ < 1, is given by 



v*(0 



1 



77 



ln(5„a ee ) + Ln(— 



(11) 



This analytical expression shows that the renormalized 
velocity in the down-spin enhances while it decreases in 



the up-spin channel. In the limit of small Eq. (11) 
is simplified and v*^(£)jv — 1 = a ee [\n(g v a ee ) =F £/2] /ir. 
All these behaviors are very familiar from the case of the 
effective mass or the effective Fermi velocity in a normal 
2DES but more significantly, the spin-polarization term 
is different than that the The discrepancy is 

due to the nature of the chiral Dirac electron behavior in 
graphene flake with having the linear dispersion relation. 
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respect to the spin direction. It would be worthwhile 
finding the asymptotic behavior of v* at some condi- 
tions. At large q, interband charge fluctuations dominate 
e(q, lo) — 1, which approaches its simple undoped system 
form. It becomes especially clear when lo is expressed in 
units of vq that the typical value of e(q, u) at large q is 
~ 1 with a non-trivial dependence on a ee . The q inte- 
grals all vary as requiring that the Dirac electrons 
model be accompanied by an ultraviolet cut-off. Since 
the crossover between intraband and interband screening 
occurs for q ~ &f, it follows that e\.£ rcs and <9 w E res have 
contributions that are analytic in a ee and vary as ln(A) 
when A is large. To leading order in ln(A) we find that 
v*/v — 1 = a ee [l - 2g v a ee g(2g v a ee )] In (A)/4 which is a 
independent and g{x) is defined in Ref.[4]. 
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FIG. 4: (Color online) Renormalized velocity of the down- 
spin, in the Dyson scheme, scaled by that of a noninteract- 
ing velocity as a function of the electron density (in units of 
10 12 cm -2 ) for different £ values at a ee = 0.5. 



In Fig. 2, we show the down-spin renormalized ve- 
locity scaled by that of a noninteracting velocity as a 
function of the coupling constant in both the dyson and 
OSA approximations which are defined by Eqs. ^ and 
respectively. Clearly, the velocity values increase 
significantly when £ approaches to unity. Despite the 
strong down-spin velocity dependence of the spin de- 
gree of freedom, the up-spin velocity becomes smoothly 

*(OSA) 

smaller with £ as it shown in Fig. 3. Notice that v a 
is always larger than v *J- Dyson \ Moreover, the £ de- 
pendence of the renormalized velocity is opposite with 
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FIG. 5: (Color online) Renormalized constant Z a for (a) the 
down-spin and (b) the up-spin as a function of the coupling 
constant a ee for cut-off value A = 100 (n = 0.36 x 10 12 cm -2 ). 



In Fig. 4 we show the down-spin renormalized veloc- 
ity as a function of the electron density (in units of 10 12 
cm 2 ). In contrast to the 2DEG, the renormalized ve- 
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locity increases by decreasing the electron density and 
indicates no Wigner crystallizatiorP^I occurs in pristine 
Dirac fermion systems^. Note that at very small n, the 
system is highly correlated and a model going beyond 
the RPA is necessary to account for increasing correla- 
tion effects at low density^. Our theoretical calculations 
show that, even at moderately low densities, the velocity 
enhancement in a supported graphene sheet can vary a 
lot in qualitatively good agreement with measured data 
in a suspended graphene sheetP^. 
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FIG. 6: (Color online) Inverse inelastic scattering lifetime 
of quasiparticle, Y1 n {k) (in units of psec -1 ) in graphene as 
a function of the on-shell energy, for different spin 

polarization values. The data in this figure refer to n = 
5 x 10 12 cm~ 2 and a ee = 0.25. 

We have also calculated the renormalization factor 
Z a {a ee , C) which is equal to the discontinuity in the mo- 
mentum distribution at k a F . The effect of £ is to make 
the values larger at large a ee compared to the case 
when £ is not included as shown in Fig. 5. The non-zero 
values of Z a , shows the Fermi liquid picture in the whole 
range of a ee and £. When n — > the Z a factor drops 
to zero logarithmically. Notice that in leading order of 
In A, the renormalization factor is independent of a and 
behaves like Z^ 1 — 1 = a ee X(2g v a ee ) In (A)/6 where X(x) 
is defined in Ref . @] . 

Finally, we compute the inelastic scattering lifetime 
of quasiparticles due to carriers-carriers interactions at 
zero temperature for different £ values. This is obtained 
through the imaginary part of the self-energy when the 
frequency evaluated at the on-shell energy 

rf„- x (k) = rf„(k) = -?3 m Ef^(k,^(k/fi)), (12) 



where Tf n (k) is the quantum level broadening of the mo- 
mentum with eigenstate |k >. ft is worthwhile to note 
that the expression of (k) is identical with a result 
obtained by the Fermi's golden rule summing the scat- 
tering rate of electron and hole contributions at wave 
vector Fig. 6 shows the behavior of the spin po- 

larization dependence of the inverse inelastic scattering 
lifetime for n = 5 x 10 12 cm~ 2 and a ee = 0.25. Imaginary 
part of the self-energy evaluated at the on-shell energy 
starts from £+(£;) = — £p, exhibits a minimum at zero 
energy and then grows up. Scattering rate in graphene is 
a smooth function which is in contrast with the conven- 
tional 2D semiconductors and 2D electron liquids because 
of the absence of both plasmon emission and interband 
processed'. We also see in Fig. 6 that the scattering 
rate is quite sensitive to the spin polarization and the in- 
elastic lifetime for minority spins is larger evidently than 
the majority spin lifetime. The ratio of the majority- to 
minority-spin lifetime is smaller than unity and related 
directly to the polarization and electron energy. 

IV. CONCLUSIONS 

In summary we have calculated the spin-polarized de- 
pendence of the quasiparticle in gaphene sheets and they 
could be strongly spin-polarization dependent and sub- 
stantially different than the usual unpolarized paramag- 
netic values. Similar to a two-dimensional paramagnetic 
diluted magnetic semiconductor electron gaJ^l, the Dirac 
electron Fermi velocity is highly spin dependent even if 
the spin polarization of the carrier population is negligi- 
bly small. Therefore, the spin-polarization dependence of 
chiral carrier transports can be observed experimentally 
nearly full spin polarization regimes. The majority-spin 
electron renormalized velocity decreases with increas- 
ing spin-polarization degree of freedom. However, the 
minority-spin electron renormalized velocity increases by 
increasing spin-polarization due to reduction in the elec- 
tron density and consequently increase in the interaction 
between electrons near the Fermi surface. We show that 
the ratio of the lifetimes of majority- to minority-spin 
electrons is smaller than unity and related to the polar- 
ization and electron energy. It has important implica- 
tions for the interpretation of many types of spin polar- 
ized experiments. Our results might be used in calculat- 
ing the effective density of states in graphene spintronic 
systems. The spin-polarized features that are the subject 



of this work may, in the future, lead to the development " MAC ALO" . 

of graphene devices incorporating interference-based spin 

filters. 
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